Strong-Coupling Theory of Counter-ions at Charged Plates 
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We present an analytical approach to the strong coupling regime of similarly and highly charged 
plates in the presence of counter-ions. The procedure is physically transparent and based on an 
exact expansion around the ground state formed by the two-dimensional Wigner crystal of counter- 
ions. The one plate problem is worked out, together with the two plates situation. Unlike previous 
approaches, the expansion is free of divergences, and is shown to be in excellent agreement with 
available data of Monte-Carlo simulations under strong Coulombic couplings. The present results 
shed light on the like-charge attraction regime. 
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The behaviour of charged particles in the vicinity of 
charged interfaces is a central yet elusive problem in the 
equilibrium statistical mechanics of Coulomb fluids, in- 
cluding colloidal science. A landmark in the field was the 
realization in the 1980s that similarly charged surfaces 
may attract each other under strong enough Coulombic 
couplings, which can be realized in practice increasing 
the valency of the counter-ions involved [H-Q ■ Notorious 
illustrations of this like-charge attraction are the forma- 
tion of DNA condensates |4j or aggregates of colloidal 
particles @. 

The weak-coupling limit is described by the Poisson- 
Boltzmann mean-field approach [6[ and by its system- 
atic improvements via the loop expansion 0-01 ■ A re- 
markable achievement of the last decade has been accom- 
plished in the opposite strong-cou plin g (SC) limit, pio- 
neered by Rouzina and Bloomfield lldl . substantiated by 
Shklovskii, Levin and collaborators Il2j , and formal- 
ized by Netz et al 13l~ll5| . An essential ingredient is that 
the layer of counter-ions close to a charged wall becomes 
two-dimensional, and in the field-theoretical method put 
forward in [14J, the leading behaviour stems from 
a single-particle theory, which produces more compact 
profiles than within mean- field theory [l6[. Next correc- 
tion orders correspond to a virial/fugacity expansion in 
inverse powers of the coupling constant S, see the defini- 
tion @ below. The method requires a renormalization of 
infrared divergences via the electroneutrality condition. 
A comparison with the Monte-Carlo (MC) simulations 
17] indicates the adequacy of the leading single-particle 
theory in the asymptotic SC limit, while capturing the 
first correction resisted the analysis. 

The establishment of an (approximative) interpolation 
between the Poisson-Boltzmann and SC regimes, based 
on the idea of a "correlation hole" , was the subject of a 
series of works Tj^- 21 [ . The specification of the correla- 
tion hole was done empirically in Refs. [l|l and self- 
consistcntly, as an optimization condition for the grand 
partition function, in [2l[ . A relevant observation in [2l[ , 
corroborated by a comparison with the MC simulations, 
was that the first correction in the SC expansion is pro- 
portional to 1/vH, and not to 1/S as suggested by the 



original SC theory. 

The aim of this Letter is to revisit the SC limit and 
establish an exact expansion which, in light of the pre- 
vious discussion, has yet to be formulated. The leading 
term of counter-ion density profiles coincides with the 
single-particle picture of the SC theory. Our expansion 
is free of infrared divergences and entails a correction in 
to the leading behaviour. Our analytical results 
are shown to be in excellent agreement with available 
MC data without adjustable parameters. The procedure 
is significantly simpler than previous works, and appears 
versatile. It will in particular be shown to yield new exact 
results in the like charge attraction regime. 

Here, we study a classical system of (equally charged) 
counter-ions in the vicinity of one or two planar walls 
bearing a uniform surface charge density, ere (e is the 
elementary charge and a > 0), the system as a whole 
being electro-neutral. The system, at thermal equilib- 
rium at the inverse temperature j3 = \/{k-QT), is im- 
mersed in a solution of dielectric constant e containing 
q-valent counter-ions, each thus having charge —qe. For 
simplicity, no image forces are present. Let us describe 
briefly the orig inal SC theory [H, G3] for the case of a 
single wall localized in the z — plane. The counter- 
ions are confined to the half-space z > 0. The rele- 
vant length scales in Gaussian units are: The Bjerrum 
length £b = (3e 2 /e, i.e. the distance at which two unit 
charges interact with thermal energy ksT, and the Gouy- 
Chapman length p = l/(27rg^Bc), i.e. the distance from 
the charged wall at which an isolated counter-ion has po- 
tential energy equal to thermal energy. All lengths r will 
be expressed in units of p, f — r/p. The counter-ion 
density profile p(z), which only depends on the distance 
from the wall z, will be considered in the rescaled form 



27rf B cr 2 



(1) 



so that the electro-neutrality condition q dzp(z) = a 
simply reads J °° dzp(z) = 1. The coupling parameter 
quantifying the strength of electrostatic correlations is 



3 = 2 7 rg 3 4c r 



(2) 
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and is large in the SC regime. According to the SC the- 
ory [ljl, [Hi , the profile of the counter-ion density can be 
formally expanded as 
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where 



h{z) 



(3) 



(4) 



The leading term po(z) comes from the single-particle 
picture of counter-ions in the linear surface-charge poten- 
tial. The MC simulations [llj indicate that the first cor- 
rection pi (z) has the expected functional form for 3 > 10, 
however, the value of the prefactor is incorrect. To be 
more particular, let us subtract the leading SC profile in 
Q and express the first correction as 



Pi(z) = 9 [p{z) ~ Po{z)} , 



(5) 



p{z) being the density profile obtained from the MC sim- 
ulations. The factor 9 can be treated as a fitting parame- 
ter which, in the original SC theory, should be equal to 3 
plus next-to-leading corrections. The numerical depen- 
dence of 9 on 3 is pictured in the inset of Fig. [1] We see 
that 6 (symbols) is much smaller than 3 (dashed line). 

Our approach is based on the fact that in the asymp- 
totic strong-coupling limit 3 — > oo, the counter- ions 
collapse on the charged surface, creating a 2D hexag- 
onal (equilateral triangular) Wigner crystal [l2| where 
every ion has 6 nearest neighbors forming an hexagon. 
Let us denote by = the position vectors 

of the vertices on this hexagonal lattice. Since there 
are just two triangles per particle, the lattice spacing 
a of the globally electro-neutral structure is given by 
q/cr = v3a 2 /2. Note that the strong-coupling limit co- 
incides with the regime in which the distance a between 
the nearest-neighbor counter-ions is much larger than the 
distance p between the counter-ions and the charged sur- 
face a = a I p, oc y/S ^> 1. In the asymptotic limit 
5 — > oo, each vertex R^ is occupied by a counter- ion i 
(i = 1, . . . , N; N —> oo). The ground-state energy of the 
counter-ion system together with the homogeneous back- 
ground charge is Eq. For 3 large but not infinite, the 
fluctuations of ions around their lattice positions start to 
play a role. 

Let us first shift one of the particles, say i = 1, from 
its lattice position Ri by a small vector SRi = (x, y, z) 
(8R\ = |<5Ri| -C a) and look for the corresponding 
change in the total energy SE = E — Eq > 0. The 
first contribution to SE comes from the interaction of 
the shifted counter-ion with the potential induced by the 
homogeneous surface charge density: SE 1 - 1 * 1 — 2-Kqe 2 az/e. 
The second contribution to 8E comes from the interac- 
tion of the particle with all other particles on the 2D 
hexagonal lattice. This contribution can be expanded as 
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FIG. 1. Comparison between the analytical first correction to 
the strong coupling profile (solid curve) and the Monte Carlo 
results of Ref. [TJ] at S = 10 3 , for a single charged wall. The 
inset compares our prediction for the rescaling factor 8 (solid 
curve given by Eq. H2J) to its Monte Carlo value reported in 
[ItJ and to the original prediction 8 = E (dashed line). 



an infinite series in x, y and z\ for our purposes, it is suf- 
ficient to consider this expansion up to harmonic terms, 
which, in the z-direction, read 
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Here, the dimensionless quantity S = X^i(^-i/ a ) 3 can 
be expressed from the general theory of lattice sums 22, 
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where ({z 7 q) = y^^ °_ n 1 / (<? + n) z is the generalized Rie- 
mann zeta function [24j and ((z) = C(z,l). Explicitly, 
5 = 11.034.... A shift of the particle simultaneously 
along all directions does not induce "mixed" harmonic 
terms of type xz or yz. The harmonic term in the (a;, y)- 
plane can be computed, and in dimensionless form, we 
have 



&5E< 
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This formula reveals a relationship between the order 
of the expansion of —f3SE in the dimensionless lengths 
x,y,z and the SC expansion in 1/vS. The linear term 
—z, which is the only one which does not vanish in the 
limit 5 — > oo, is the leading term. It corresponds to the 
single-particle picture, in close analogy with the previ- 
ous SC theory. The harmonic terms turn out to be of 
the SC order (3(qe) 2 p? /a 3 oc l/\/3 and likewise, terms of 
the pth order in the variables x, y, z are of the SC order 
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(3(qe) 2 p p /a p+1 oc l/S^ 1 )/ 2 . This scheme constitutes a 
systematic basis for SC expansions. 

The generalization of the above formalism to all par- 
ticles is straightforward. We shift every particle i — 
1,2, ... ,N from its lattice position by a small vec- 
tor SRi = (xi,yi, Zi). In what follows however, we shall 
be interested in the counter-ion density profile which only 
depends on the z coordinate. Thus, when expanding in 
statistical averages the Gibbs weight exp(— (35E) in pow- 
ers of 1/ we can restrict ourselves to the z-harmonic 
part. The corresponding change in the total energy 5E 
is given by a counterpart of ([8]), 
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The next simplification comes from the fact that parti- 
cles are identical, exposed to the same potential induced 
by the surface charge, so that a summation over parti- 
cle coordinates can be represented by just one auxiliary 
coordinate. For the density particle profile, defined by 
P( r ) = (Z)i=i <K r - r i)) = N(S(r- ri)), we get explicitly 



p{z) = Ce~ s 



dz'e 



1 + 
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where C is determined by the normalization condition 
/ p = 1. Simple algebra gives 
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Comparing this result with the previous one ([3]), Q ob- 
tained in the original SC theory [l3, 14 1, we see that the 
leading terms coincide, while the first corrections have 
the same functional dependence in space but different 
prefactors. The result (jlip can be re-expressed in terms 
of the 6- factor, introduced in the relation ([5]), as follows 



( 47r ) 3 / 2 1 = 



(12) 



This formula, in excellent agreement with MC data, dif- 
fers substantially from the previous SC estimate 9 = S, 
see Fig. [TJ 

The method can be readily applied to the case of two 
parallel walls, each having the same charge density ere, 
located at distance d from one another. The electric field 
between the walls is equal to now. At T — 0, the clas- 
sical system is defined by the dimensionless separation 
r] = d\J a jq = d/y/2nE. A complication comes from the 
fact that counter-ions form, on the opposite surfaces, a 
bilayer Wigner crystal, the structure of which depends 
on rj Five different structures are energetically 

favored for various regions of rj. At the smallest separa- 
tions when rj — > 0, the single hexagonal lattice (structure 



I, see Fig. [2]- left), with rows distributed consecutively 
between the two surfaces, is formed. At large separa- 
tions r/ > 0.732, the energetically favored geometry is 
composed of two staggered hexagonal lattices (structure 
V), one for each plate. We shall document our SC ap- 
proach on structure I, relevant to obtain the large S be- 
haviour. Due to global neutrality, the lattice spacing b 
of the single (bilayer) hexagonal structure is given by 
q/(2a) = V3b 2 /2. The SC regime is identified with the 
condition d <C b, or equivalently 5> d 2 . 
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FIG. 2. Left: Structure I of counter-ions on two parallel 
charged plates (see text). Right: Phase diagram following 
from the equation of state (|17p : the solid curve shows the 
points where P = 0. The filled squares are those MC results 
from Ref [l7| ] with S > 20, while the dashed line is for the 
points where dP/dd = 0, which defines d m ax. 



The two walls are located at positions z = and z = d. 
The position vector of the particle localized on the 
shared hexagonal Wigner lattice will be denoted as R^ 
if it belongs to the wall at z — (say filled symbols of 
the left panel of Fig. [2]) and as Ii\ if it belongs to the 
wall at z — d (open symbols in Fig. [5]). Let us shift 
the particle i — 1 localized on the z — wall by a small 
vector JR^ = (x, y, z) and look for the energy change 
5E from the ground state. Since the potential induced 
by the surface charge on the walls is constant between 
the walls, the corresponding SE^ = 0. The harmonic 
term in the z-direction reads 



eSE^ = 
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Using the exact values of the partial hexagonal sums 



E^i b/R. 



= 75/12, SE { Z 



? (0) 1 = 55/12, £ 4 [b/Rf 
turns out to be positive, as it should. The harmonic term 
in the {x, y)-plane can again be computed but proves im- 
material for the sake of our purposes. When all par- 
ticles are shifted from their lattice positions {R;} by 
{<5Ri = (x{,yi, z^}, the total energy change is given, as 
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far as the z-dependent contribution is concerned, by 

(t 2 
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Expanding exp(— in 1/yH and enforcing electro- 
neutrality, the density profile p(z) is obtained in the form 
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where 9 
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This 9 differs from the single-plate one (IT21 by the factor 
l/y/2 due to the different hexagonal lattice spacings a 
and b. The functional form of (fT5|) coincides with that 
of Moreira and Netz TH, 14 1. For (not yet asymptotic) 
H = 100, the previous SC result 9 = S is far away from 
the MC estimate 9 ~ 11.2 [17], while our formula (|16|) 
gives ~ 12.5. 

Applying the contact-value theorem to the density pro- 
file (|T5j) . the pressure P between the plates is given by 



P 



P 



2 d ^ d 2 

J + aS + ?. 
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An analogous result was obtained within the approximate 
approach of Ref. 21], with the underestimated ratio 
0/y/E = 3^3/2 = 0.866 .... We recall that our estimate 
of 9 is valid only in the structure I region < 77 < 0.006; 
increasing 77, other energetically favored bilayer Wigner 
structures have to be considered as a starting point. Eq. 
(fTTf provides insight into the like charge attraction phe- 
nomenon. The attractive (P < 0) and repulsive (P > 0) 
regimes are shown in Fig. [2] (right panel). Although our 
results hold for d -C S 1 / 2 and for large S, the shape of 
the phases boundaries where P — (solid curve) shows 
striking similarity with its counterpart obtained numer- 
ically [171 Eil ] (we note for instance that the terminal 
point shown by the filled circle in Fig [2] is located at 
d = 4, a value close to that which can be extracted 
from [It], [3]). While the upper branch of the attrac- 
tion/repulsion boundary is such that d/y/E. is of order 
unity and hence lies at the limit of validity of our ex- 
pansion, we predict the maximum attraction to be ob- 
tained for (i m ax = v60 cx S 1 / 4 , as follows from enforcing 
dP/dd = 0. Since d ma ^/\/E cx S~ 1//4 — > 0, we can con- 
sider the latter prediction, shown by the dashed line in 
Fig. [2] as asymptotically exact; we note that it is fully 
corroborated by the scaling laws reported in (l^ |. 

In conclusion, the present exact approach to the strong 
coupling regime shows that while the leading order re- 
sults at large S (for one or two plates) can be obtained 



by a single counter-ion theory, the next terms actually re- 
flect the complete ground state structure (N counter-ion 
property). This explains the failure of virial-like expan- 
sions. We have shown how such shortcomings can be 
circumvented within a physically transparent procedure, 
and obtained analytical results in remarkable agreement 
with Monte Carlo data. In practice, for a highly charged 
interface in water at room T, one has ai^ ~ 1, so that S 
takes values close to 50, 170 and 400 for respectively di-, 
tri-, and tetra-valent counter-ions, as often found in bi- 
ology (spermine). Although asymptotic, our predictions 
turn out to be reliable for such couplings. A generaliza- 
tion of the approach to dielectric inhomogeneities [28| . 
systems with salt or asymmetric [2^ |. and curved sur- 
faces (3(|, offer interesting problems for more detailed 
studies in the future. The formulation is also convenient 
for a quantum-mechanical generalization. 
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